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Abstract. We study some functional inequalities satisfied by the distribution of the 
solution of a stochastic differential equation driven by fractional Brownian motions. Such 
functional inequalities are obtained through new integration by parts formulas on the 
path space of a fractional Brownian motion. 
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1. Introduction 
Let (Xf )t>o be the solution of a stochastic differential equation 



where {Bt)t>o is a ri-dimensional fractional Brownian motion with Hurst parameter H > ^. 
Under ellipticity assumptions and classical boundedness conditions (see ^ and [H]), the 
random variable X^, t > 0, admits a smooth density with respect to the Lebesgue measure 
of MJ^ and the functional operator 

pj{x) = E{f{xn) 

is regularizing in the sense that it transforms a bounded Borel function / into a smooth 
function Ptf for t > 0. In this note we aim to quantify precisely this regularization 
property and prove that, under the above assumptions, bounds of the type: 

\V,, ■ ■ ■ V,,Ptf{x)\ < Ci,...,,{t,x)\\f\U t > 0, X G M", 
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are satisfied. We are moreover able to get an explicit blow up rate when t — )■ 0: For a 
fixed X e M", when t — )■ 0, 

Our strategy to prove such bounds is the following. If / is a C°° bounded function on 
M", we first prove (see Lemma [4. ip that the following commutation holds 

V,Ptf{x) = E {^a',{t,x)Vtf{Xt)^ , 

where the a{t, x)'s solve an explicit system of stochastic differential equations. Then, using 
an integration by parts formula in the path space of the underlying fractional Brownian 
motion (see Theorem Theorem 13.61) we may rewrite the expectation of the right hand 
side of the above inequality as E ($i(t, x)/(Xf )) where $j(t,x) is shown to be bounded 
in L^, 1 < p < +00 with a blow up rate that may be controlled when t 0. It yields a 
bounds on \ViPtf{x)\. Bounds on higher order derivatives are obtained in a similar way, 
by iterating the procedure just described. Let us mention here that the bounds we obtain 
depends on bounds for the inverse of the Malliavin matrix of Xf . As of today, to 
the knowledge of the authors such bounds have not yet been obtained in the rough case 
H < ^. The extension of our results to the case H < ^ is thus not straightforward. 

We close the paper by an interesting geometric situation where we may prove an optimal 
and global gradient bound with a constant that is independent from the starting point 
X. In the situation where the equation is driven by a Brownian motion such global 
gradient bound is usually related to lower bounds on the Ricci curvature of the Riemannian 
geometry given by the vector fields V^'s, which makes interesting the fact that the bound 
also holds with fractional Brownian motions. 

2. Stochastic differential equations driven by fractional Brownian 

MOTIONS 

We consider the Wiener space of continuous paths: 

W" = (C([0,l],M"),(i30o<t<i,P) 

where: 

(1) C([0, 1],M") is the space of continuous functions [0, 1] — )■ M"; 

(2) Wt)t>o is the coordinate process defined by f3t{f) = f (t), f G C([0, 1],M"); 

(3) P is the Wiener measure; 

(4) {Bt)o<t<i is the (P-completed) natural filtration of (/3t)o<t<i- 

A n-dimensional fractional Brownian motion with Hurst parameter H G (0, 1) is a Gauss- 
ian process 

Bt = iBl...,B^), t>0, 
where B^, . . . , i?" are n independent centered Gaussian processes with covariance function 

It can be shown that such a process admits a continuous version whose paths are Holder 
7 continuous, 7 < if. Throughout this paper, we will always consider the 'regular' case. 
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H > 1/2. In this case the fractional Brownian motion can be constructed on the Wiener 
space by a Volterra type representation (see [3]). Namely, under the Wiener measure, the 
process 

(2.1) Bt= f KH{t,s)dP,,t>0 

Jo 

is a fractional Brownian motion with Hurst parameter H, where 

1 /"* 3 1 

KH{t,s) = chs^~^ {u - s)^~^u^~^du , t>s. 
and ch is a suitable constant. 

Denote by £ the set of step functions on [0, 1]. Let "H be the Hilbert space defined as the 
closure of £ with respect to the scalar product 

(l[o,i], l[o,s])w = Rnit, s). 
The isometry from to -^^^([0, 1]) is given by 

{K*^^){s) = l\{t)^{t,s)dt. 

Moreover, for any ip G L^([0, 1]) we have 

[\{s)dB,= f\Kl^){s)d(5,. 
Jo Jo 

Let us consider for x G M" the solution {X^)t>o of the stochastic differential equation: 



(2.2) 



i=i -^0 



where the V^'s are C°° bounded vector fields in R". Existence and uniqueness of solutions 
for such equations have widely been studied and are known to hold in this framework (see 
for instance [10]). Moreover, the following bounds were proved by Hu and Nualart as an 
application of fractional calculus methods. 

Lemma 2.1. (Hu-Nualart, Consider the stochastic differential equation Ii2.2\) . If the 

derivatives of Vi 's are bounded and Holder continuous of order A > 1/H — 1, then 



E ( sup \Xt\A 

\0<t<T J 



< OO 



for all p > 2. If furthermore Vi 's are bounded and E(exp(A|Xo|'')) < oo for any A > and 
q < 2H , then 



i exp A ( 



E ( exp A ( sup \Xt\'' ] ] < oo 

.0<t<T 



for any A > and q < 2H . 



Throughout our discussion, we assume that the following assumption is in force: 
Hypothesis 2.2. 

(1) Vi{x) 's are bounded smooth vector fields on R" with bounded derivatives at any order. 

(2) For every x G R", (Vi(x), ■ ■ ■ , Vn{x)) is a basis o/R". 
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Therefore, in this framework, we can find functions cot such that 



(2.3) [^.,^.] = E4^^ 



k=\ 

where the w^^'s are bounded smooth functions on with bounded derivatives at any 
order. 



3. Integration by parts formulas 
We first introduce notations and basic relations for the purpose of our discussion. Consider 

dx 



the diffeomorphism x) = Xf : M" — )■ M". Denote by Jt = the Jacobian of ■). 



It is standard (see [IT] for details) that 

(3.4) dit = ^dVi{X't)itdBl with Jo = I, 

i=l 

and 

n 

(3.5) c^Jj-^ = - ^:i^^^Vi{X't)dBl with Jq^ = I. 

i=l 

For any vector field W on R", we have that 

{^tM){Xf) = itW{x), and {^t-^W){x) = Jr'W^(Xn 

Here ^t* is the push-forward operator with respect to the diffeomorphism $(/!:, x) : — )■ 
M". Introduce the non-degenerate n x n matrix value process 

(3.6) a{t,x) = {ai{t,x))l^^, 
by 

n 

($i,V,)(Xn = Jtmx)) = J2c^lit,^)Vk{Xn ^ = 1,2, ...,n. 

fc=i 

Note that a{t, x) is non-degenerate since we assume l^'s form a basis at each point x G M". 
Denote by 

(3.7) I3{t,x) = a^^{t,x). 
Clearly we have 

n 

(3.8) ($,;V,(X,^))(x) = JrV,(Xn(x) = ^/3^(t,x)\4(x) ^ = 1,2,.. .,n. 

k=l 

Lemma 3.1. Let a(t,x) and f3{t,x) be as above, we have 

n 

(3.9) dai{t,x) = -J2 (^kit,xWikiXt)dBl, with aj(0,x) = S}; 

k,l=l 

and 

n 

(3.10) d(3}{t,x) = J2 0J^.{XnPf{t,x)dBl, with /3j(0,x) = 6]. 

kl=l 
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Proof. The initial values are apparent by the definition of a and /3. We show how to 
derive equation fl3.10p . Once the equation for f3(t,x) is obtained, it is standard to obtain 
that a(t,x) = f3~^(t,x). 

Consider the n x n matrix V = (Vi, V2, Vn) = (Vj) obtained from the vector fields V. 
Let W be the inverse matrix of V. It is not hard to see we have 



[X . 



/3i(t,x) = ^l^^-(x)(J-V,(Xf))^ 

k=l 

By the equation for Xf , relation fl2.3p . equation (13. Sp . and Ito's formula, we obtain 

n 

k=l 

n 

= Y.^i^(^t){3t'Vl{XnKx)dB^. 



Hence 



□ 



n. 



dp%x)=Y,ujl{X^)fi%x)dBl. 

k,l=l 

This completes our proof. 
Define now hi{t, x) : [0,1] x W ^ K hj 

(3.11) h,it,x) = {(3'[is,x)I[o,t]{s))k=i,...,n, I = 1 
Introduce Mij(t,x) given by 

(3.12) M,,,(t,x) = :^{hi{t,x),hj{t,x))n. 
For each t G [0, 1], consider the semi- norms 

11,11 l/(^)-/(^)l 
Wflkt ■= sup — -— . 

0<v<u<t [U - Vp 

The semi- norm ||/||7,i will simply be denoted by ||/||7. 
We have the following two important estimates. 

Lemma 3.2. Let a{t,x), P{t,x) and hi{t,x) be as above. We have: 

(1) For any multi-index v, integers k,p> 1, there exists a constant Ck^p{x) > depending 

on k,p and x such that for all x G M" 



sup 

0<t<l 



dx" 



a{t, x] 



< Ck,p{x 



k,p 



sup 

0<t<l 



dx'^ 



P{t,x) 



< Ck,p{x). 



k,p 



(2) For all integers k,p > 1, 6hi{t,x) G D*^'^. Moreover, there exists a constant Ck,p{x) 
depending on k,p and x such that 

\\6hi{t,x)\\k,p<Ck,p{x)t'', te [0,1]. 

In the above S is the adjoint operator o/D. 
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Proof. The result in (1) follows from equation (13.91) . (13.101) and Lemma 12.11 In what 
follows, we show (2). Note that we have (c.f. Nualart|;9j) 



6hi(t, x) 



Jo Jo Jo 



h(t, x)^\u — V 



2H-2 



dudv 



f3i{u,x)dBu - an 



t rt 



^0 



DuPiiv, x)\u — v\ dudv 



Here an = H{2H — 1). From the above representation of 6hi and the result in (1), it 
follows immediately that 6hi{t,x) G D'^'^ for all integers k,p> 1. To show 

\\Shi{t,x)\\k,p < Ck,p{x)t^ for all t G [0, 1], 

it suffices to prove 



Pi{u,x)dBu 



<C{x)t^ tG[0,l]. 



Here C{x) is a random variable in L'^{F). Indeed, by standard estimate, we have 



(A(n,x) - A(o,x))d5„ 







<Cm-,x)\\r,t\\B\U,t, tG [0,1]. 



In the above | < r, 7 < if and r + 7 > 1, and C > is a constant only depending on 7. 
Therefore 



Pi{u,x)dBu 



< C\\P{-,x)\\r^t\\B\\,^tt^ + \P{0,x)\\Bt\, t G [0, 1]. 



Together with the fact that for any t < H, there exists a random variable Gr{x) in L^(P) 
for all p > 1 such that 

\l3{t,x) - I3{s,x)\ < Grix)\t - sT, 

the proof is now completed. 



□ 



Lemma 3.3. Let M{t,x) = {Mij{t,x)) be given in li3.1S\) . We have for all p > 1, 

sup E [det(M(t,x))^*'] < 00. 

t€[0,l] 

Proof. Denote the Malliavin matrix of Xf by r{t,x). By definition 



Tij{t,x) = {-D,Xi,-D,Xi)n = an 



t rt 



Jo 



T>uXi'D^Xl\u - V 



2H~-2 



dudv. 



It can be shown that for all p > 1 (cf. Baudoin-Hairer [3j, Hu-Nualart [7j and Nualart- 
Saussereau [llj). 



(3.13) 

Introduce 7 by 



sup E ( det 

te[o,i] 



T{t,x] 

t2H 



< 00. 



%,j{t,x) = an 



t rt " 



JO 



J2 {JZ'Vk{Xu)y {3-'V,iX.,)Y \u - v^-'dudv. 



k=l 



Since Df'Xi = JtJ^ ^Vk{Xs), we obtain 



(3.14) 



T{t,x) = 3a(t,x)3i 
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Recall 

Mij{t,x) = ^{hi{t,x),hj{t,x))n, 

where 

hi{t,x) = {(3f{s,x)Iio^t]{s))k=i,...,n, i = l,...,n. 
By (ESD and (KIM , we have 

(3.15) Vix)Mit,x)V{xf = :^l{t,x) = J-i£M(j-i)T. 

Finally, by equation (13 ■4p . Lemma [2. II and estimate (I3.13P we have for all p > 1 

sup E [det(Mij)"^] 

< oo, 

te[o,i] 

which is the desired result. □ 
The following definition is inspired by Kusuoka |8]. 

Definition 3.4. Let H be a separable real Hilbert space and r G M &e any real number. 
Introduce ICr{H) the set of mappings $(t, x) : (0, 1] x R"- — )■ I}°°[H) satisfying: 

(1) ^(t, x) is smooth in x and ^^{t,x) is continues in {t,x) G (0, 1] x M" with probability 
one for any multi-index v; 

(2) For any n,p > 1 we have 

sup t 



0<t<l 

We denote /C,(M) by /C^. 



d^x 



{t,x] 



< oo. 



Lemma 3.5. With probability one, we have 

(1) a(t,x), f3(t,x) G /Co; 

(2) 6hi{t,x) G /Ci; 

(3) Let [M^^) be the inverse matrix of (Mij). Then M^j- G /Co for all i,j = 1, n. 

Proof. The first two statements are immediate consequences of Lemma 13. 2[ The third 
statement follows by writing = estimates in Lemma 13.21 (1) and Lemma 

1331 ' □ 

Now we can state one of our main results in this note. 

Theorem 3.6. Let f be any C°° bounded function and $(t, x) : Q ^ ICr we have 

E mt, x)vj{xn) = E ((r^^$(t, x))f{xn) , 

where Ty^$(t, x) is an element in ICr-i with probability one. 

Proof. This is primarily integration by parts together with the estimates obtained before. 
First note 

Di/(X,) = (V/(X,),D^,X,) 

= (v/(x,),JtJ;V,-(x,)) 

n 
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Hence 

(3.16) VJ{X,) = -i^ f^/3i(t)M-i(D/(X,),/.Ki))«- 

^ 3,1=1 

Therefore, we have 

E($(t,x)V,/(Xi)) 
1 " 

=:^Y.^{{'Df{x^)Mt.^)mM,,\t)h,{t))H) 

kl=l 



^2H 

k,l=l 



k,l=l 



By Lemma [3.51 the first term in the brackets above is in /C^-i and the second term is in 
Kr- Finally, denote 



T*$(t,x)= ^ 



kl=l 



It is clear that Ty.$(t,x) G /Cr-i- The proof is completed. □ 

4. Gradient Bounds 

With the integration by parts formula of Theorem 13.61 in hands we can now prove our 
gradient bounds. We start with the following basic commutation formula: 

Lemma 4.1. For i = 1,2, n, we have the commmutation 

\k=l 

where the a{t,x) solve the system of stochastic differential equations ^3.9^) . 
Proof. For any C^-vector field W on M" we have 

l^P,/(x)=E(((W)/)(Xf)). 
The remainder of the proof is then clear from the computations in the previous section. □ 
Finally we have the following gradient bounds. 

Theorem 4.2. Let p > 1. For ii, ...,ik G {1, n}, and x G, M", we have 

\Vi,...Vi,Ptf{x)\ < c{t,x){pjp{x)y^ t G [0,1] 

with C(t,x) = O (jwk) when t — )• 0. 
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Proof. By Theorem 13.61 and Lemma for each k >1 there exists a ^^^''\t^x) G /C_fe 
such that 

V,,...\/,,Pi/(x) =E($(-^')(t,x)/(Xj). 
Now an apphcation of Holder's inequahty gives us the desired result. □ 

Remark 4.3. Here let us emphasize a simple but important consequence of the above 
theorem that, suppose f is uniformly bounded, then 

\V,,...V,,Ptf{x)\<C{t,x)\\f\\^ tG [0,1] 

where C(t,x) = O {jet;) as t ^ 0. 

Another direct corollary of Theorem 14.21 is the following inverse Poincare inequality. 
Corollary 4.4. For ii, ...,ik G {1, n}, and x G M", 

\Vi,...V,^PJ{x)f < C{t,x){Pj'{x) - {Ptf)\x)) t G [0,1] 
with C(t,x) = O (t^^t) when t 0. 

Proof. By Theorem 14.21 for any constant C G M we have 

\V,,...V,^PJix)\' = \V,,...V,Mf - C)ix)\' < C{t,x){Pt{f - Cf{x)) t G [0, 1] 

with C(t, x) = O {-jpTk) when t — )■ 0. Now minimizing C G M gives us the desired 
result. □ 

Remark 4.5. For each smooth function / : — )■ M, denote 

n 

r(/) = 

i=l 

We also have, for ii, i^. G {1, n}, and x G M", 

\V,,...V,,PJ{xf <C{t,x)PtT{f){x), tG [0,1] 
with C{t,x) = O ( ^2g(fc-i) ) when t — > 0. Indeed, by Theorem \3.6\ and Lemma 4-P we know 
that for each k >1, there exists ''\t,x) G /Ci_fc,j = 1,2, ...,n such that 

V,,...V,,Ptfix) = E ($(^-')(t,x)(V,/)(X,)) . 
The sequel of the argument is then clear. 

5. A GLOBAL GRADIENT BOUND 

Throughout our discussion in this section, we show that under some additional conditions 
on the vector fields Vi, ...,Vn, we are able to obtain 

uniformly in x, where we denoted as above 

n 
i=l 

For this purpose, we need the following additional structure equation imposed on vector 
fields V^„...,Ki. 
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Hypothesis 5.1. In addition to Hypothesis \2.S\ we assume the smooth and bounded 
functions cjf- satisfy: 



k j 



I < i, j, k < d. 



Interestingly, such an assumption already appeared in a previous work of the authors |1] 
where they proved an asymptotic expansion of the density of Xt when t — )■ 0. 

Remark 5.2. In the case of a stochastic differential equation driven by a Brownian 
motion, the functional operator Pt is a diffusion semigroup with infinitesimal generator 
L = \ (^"=1 '^i)- gradient subcommutation 

v/f(iV) <P,(v^), 

is then known to be equivalent to the fact that the Ricci curvature of the Riemannian 
geometry given by the vector fields Vi's is non negative (see for instance 

The foUowing approximation resuh, which can be found for instance in [6j, will also be 
used in the sequel: 

Proposition 5.3. For m > 1, let = {-BJ"; t G [0,1]} be the sequence of linear 
interpolations of B along the dyadic subdivision of [0, 1] of mesh m; that is if t™- = i2~^ 
fori = 0, ...,2'^; then fort G (tr^^ml' 

■j^m _ _| 1 (Bfm — Bf-mY 

Consider X"^ the solution to equation ^2.2^ restricted to [0, 1], where B has been replaced 
by B™' . Then almost surely, for any 7 < and t G [0, 1] the following holds true: 



(5.17) 



lim IIX^'-X'"' 







Theorem 5.4. Recall the definition ofa(t,x) in jiS. 6\) and 

n 

(5.18) rfa}(t,x) = - ^ ai{t,x)ui{X^)dBl, with a5(0,x) = S]. 



kl=l 



Under Assumption \5 . 1\ uniformly in t and x, we have 



(5.19) 



almost surely. 



X] 



< 1; 



and 



X) 



i=l 



< 1, 



Proof. Let us thus consider and a'^{t,x) the solution of (12.21) and fl5.18p where B is 
replaced by B"^, that is 



dxr = J2v,{xr)dBT'\ 

i=l 

n 

da"'{t,x) = -^a'"(t,x)wfc(Xr)c/F 



m,k 
s 1 



k=l 
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with Xq = X and a(0,a;) = I. Here = {^kj)- order to show that the process a{t,x) 
is uniformly bounded, by applying Proposition 15.31 to the couple (X, a), it is sufficient to 
prove our uniform bounds on a;™(t,a;), uniformly in m. In the sequel set 



k,m 



r)k,m r)k,m 

''n '■n-l 

j-m j-m 



for 1< n < 2"^ and 1< < n. 



Then, for t G [C-i,C), we have 

n 

Therefore, for t G we obtain 



k=l 



a^{t,x) 



n— 1 n n— 1 



Proceeding inductively, we end up with the following identity, valid for t G [t™-!,^™) and 



n = 0,...,2'": 



(5.20) a"'{t,x) 



By Assumption 15.11 each Uk is a skew-symmetric matrix, expression (15.201) gives us 



«j (t, 



< 1; and 



X 



i=l 



< 1. 



This is our claimed uniform bound on a"^{t, x), from which the end of our proof is easily 
deduced. 

□ 



As a direct consequence of Lemma 14.11 and Theorem 15.41 we have the main result of this 
section. 



Theorem 5.5. Under A s sumption \ 5. 1[ we have uniformly in x 

Vnpj)<Pt{Vm). 



Proof. By applying Lemma 14. H Cauchy-Schwarz inequality and then Theorem 15. 4[ we 
have for any vector a = (oj) G M" 



J2aiViPtf{x) =eIJ2 a,al{t,x)Vkf{X, 



i=l 



Ki,k=l 



< E 



n In 



2\ 2 



Y.(^kf{xnf 



k=l \i=l 



.fc=l 



< \\a\\Elj2Vkf{X^) 
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By choosing 



{ViPtf){x) 



we obtain 



Vmf) 



The proof is completed. □ 

Remcirk 5.6. Since Pt comes from probability measure, we observe from Jensen inequality 
that 

Vmr)<PtiVm) 

implies 

nPtf) < PtiTf). 
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